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ABSTRACT 

To study phenomena of plasmas around rotating black holes, we have derived 
a set of 3+1 formalism of generalized general relativistic magnetohydrodynamic 
(GRMHD) equations. Especially, we investigated general relativistic phenomena 
with respect to the Ohm's law. We confirmed the electromotive force due to 
the gravitation, centrifugal force, and frame-dragging effect in plasmas near the 
black holes. These effects are significant only in the local small-scale phenomena 
compared to the scale of astrophysical objects. We discuss the possibility of 
magnetic reconnection, which is triggered by one of these effects in a small- 
scale region and influences the plasmas globally. We clarify the conditions of 
applicability of the generalized GRMHD, standard resistive GRMHD, and ideal 
GRMHD for plasmas in black hole magnetospheres. 

Subject headings: plasmas, general relativity, methods: analytical, galaxies: ac- 
tive, galaxies: jets, galaxies: magnetic fields, galaxies: nuclei 



1. Introduction 

Numerical simulations of general relativistic magnetohydrodynamics (GRMHD) have 
revealed a number of interesting and important physics of plasmas in black hole magneto- 



quasars (uQSOs), and gamma-ray bursts (GRBs) ( 


K^oide et al. 


1998. 


1999, 


2000; 


Koide 


2004; 


Koide et al. 


2006: 


Gammie et al. 


2003; 


McKinnev 


2006). All of these GRMHD simulations 



were performed within an assumption of zero electric resistivity (ideal GRMHD). An order es- 
timation of the global plasma variables with respect to accretion disks ar ound almost all kinds 
of black holes suggested validity of the ideal GRMHD (jMcKinneyll2004l ) . On the other hand. 
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all long-term GRMHD simulations of jet formation in black hole magnetospheres showed 
artificial appearance of magnetic islands, which are caused through magnetic reconnections 
due to numerical resistivity. In spite of the numerical inconsistency, these numerical results 
clearly suggested spontaneous formation of anti-parallel magnetic configuration, where mag- 
netic reconnection is caused easily in the black hole magnetospheres. The magnetic reconnec- 
tion would change the global magnetic configuration drastically and influence the global dy- 
namics of plasmas around the black holes. Thus, calculations including resistivity, the cause 
of magnetic reconnection, are required. In this aspect, special relativistic magnetohydrody- 
namics (sRMHD) with electric resi stivity (resistive sRMH D) has been utilized to mimic the 
relativistic magnetic reconne ction (IWatanabe et al.ll2006f). Also, resistive GRMHD has be en 
discussed by several autho rs ( Bekenstein k. Oron IQTsI : Khanna k. Camenzind 1994 . 1996al lbl: 
Kudoh &: Kaburakilll996l ). and applied to accretion disks around Kerr (i.e., rotating) black 
holes. In both resistive sRMHD and GRMHD, the authors have used the standard Ohm's 
law. In spite of the mathematical consistency of the resistive sRMHD and GRMHD with the 
standard Ohm's law (standard sRMHD/GRMHD), we should use the results of these calcu- 
lations carefully, because causality is bro ken and artif icial wave instability is caused because 
of the usage of the standard Ohm's law (iKoidd 120081 ). To guarantee causality with electric 
resistivity, we have to use generalize d sRMHD or GRMHD including the generalized rela- 
tivistic Ohm's law (IKoidd l2008l . l2009l ). The generalized GRMHD equations were introduced 
on the basis o f the two-fluid approximation of plasma in the Kerr metric by the pioneer, 
Khannal (119981 ). More generalized equations from the ge neral r e lativis tic Vlasov-Boltzmann 
equation in time- varying space-time were formulated by iMeierl (120041 ). With respect to the 
generalized sRMHD equations derived from relativistic two-fluid equations, it was proved 
that c ausality is sa ti sfied fo r the p air plasma whose plasma parameter is much greater than 
unity (IKoidd I2OO8I ). iKoidd (120091 ) extended these generahzed sRMHD equations of a pair 
plasma to those of any two-component plasmas including not o nly the pair plasma but also 



the electron-ion plasma. These generalized sRMHD equations of lKoidd (12008 
special relativistic basic phenomena of plasmas. 



20091 ) revealed 



In this paper, we extend the generalized sRMHD equations of iKoidd (l2008l . 120091 ) to 
general relativistic version to investigate the distinctive phenomena of plasmas in the black 
hole m agnetospheres. Com paring t he ge neralized GRMHD equations suggested by iKoide 
( 120091 ) with those derived by iMeien (120041 ) , we found that we should not use the assumption 
of the infinitely small difference of the variables with respect to the enthalpy introduced 
by iKoidd (120091 ). We also found that the condition of 4- velocity, which is a null vector, 
is automatically satisfied when we use the appropriate definitions of the mass density and 
4-velocity for the one-fiuid approximation. Now, the genera l ized G RMHD equations are 
derived without the sever restrictions as those used in iKoidd ( 120091 ) . This means that the 
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generalized GRMHD equations presented here are identical to the general relat ivistic two- 



fluid equations mathematically. Concerning the comparison with iMeierl (120041). w e foun d 
good correspondence between our GRMHD equations and those derived by iMeierl (12004 ). 
while they are not identical (see Section H]). To clarify and evaluate the distinctive nature 
of the plasmas around black holes, we derive a 3+1 formalism of the generalized GRMHD 
equations in a fixed space-time around rotating black holes. Especially, we concentrate on 
the distinctive properties suggested by the generalized Ohm's law of plasmas around rotating 
bl ack holes . We a lso found that this 3+1 formalism corresponds to the equations derived 



by iKhannal (Il998l ) excellently when we consider a cold plasma, where the pressure is much 
smaller than the rest mass energy density. The 3+1 formalism of the generalized GRMHD 
equations will be useful when we perform numerical simulations of the plasmas around black 
holes, including the resistive and Hall effects within causality. 

In Section [2], we derive the generalized GRMHD equations and their 3+1 formalism 
based on the general relativistic two-fluid equations. We clarify the distinctive phenomena of 
the plasmas around rotating black holes and the conditions for applicability of the generalized 
GRMHD, standard resistive GRMHD, and ideal GRMHD equations in Section [3l The last 
section presents discussions. 



2. Generalized GRMHD equations 

2.1. Covariant form 

We derive generalized GRMHD equations based on the general relativistic two-fluid 
equations. For simplicity, we assumed that the plasma is composed of two fluids, where one 
fluid consists of positively charged particles with mass m+ and electric charge e, and the other 
fluid consists of negatively charged particles with mass m_ and electr ic charge — e. Unlike 



the discussion of the generahzed GRMHD equations in iKoidd (120091 ) . we use the general 
relativistic two- fluid equations without simplification, i.e., without the conditions of non- 
relativistic relative velocity of the two fluids of the plasmas, of non-relativistic pressure, and of 
negligible difference of a certain normalized enthalpy of the two fluids. We take no account of 
radiation cooling effect, plasma viscosity, and self-gravity in order to study the fundamentals 
of interaction between magnetic fields and resistive plasmas around the rotating black holes. 
We also assumed that the plasmas are heated only by Ohmic heating and disregarded nuclear 
reactions, pair creation, and annihilation. We neglected quantum effects of the elemental 
processes in the plasmas. The space-time, {x^ , , x"^ , x^) = {t,x^,x'^,x^), is characterized by 
a metric gf^,^, where a line element is given by ds"^ = g^^dx'^dx'^ . Here, we use units in which 
the speed of light, the dielectric constant, and the magnetic permeability in vacuum all are 
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unity: c=l, eo = l, /io = l- The relativistic equations of the two fluids and the Maxwell 
equations are 

V.(n±f/^) = 0, (1) 
V,(/i±f/^f/^) = - V>± ± en±(?'^'^[/^F., ± i?^ (2) 
V, = 0, (3) 

where variables with subscripts, plus (+) and minus (-), are those of the fluid of positively 
charged particles and of the fluid of negative particles, respectively, n± is the proper particle 
number density, p± is the proper pressure, h± is the relativistic enthalpy densitjo, is the 
4-velocity, is the covariant derivative, F^^, = V^A^, — V^A^ is the electromagnetic field 
tensor (A^ is the 4- vector potential), *F^'^ is the dual tensor of F^^,, is the frictional 4- 
force density between the two fluids, and J'^ is the 4-current density. We will often write a set 
of the spatial components of the 4-vector using a bold italic font, e.g., U± = (U]_, f/|, f/|), 
J = {J\ J2, .P), R= {R\R^,R^). We further define the Lorentz factor 7± = U^, the 3- 
velocity = U^/'y±, the electric field Ei = F^\ the magnetic flux density J2k=i ^ijkBk = 
[eijk is the Levi-Civita tensor), and the electric charge density pe = J^- Here, the alphabetic 
index (z, j, k) runs from 1 to 3. 

To derive one-fluid equations of the plasma, we define the average and difference vari- 
ables as follows: 

p = -|- m_n_7^_, (5) 

n = ^, (6) 
m 

P = P+ + P-, (7) 

Ap = P+-P-, (8) 

W = -(m+n+U^ + m^n^U'l), (9) 

JA* = e(n+t/^ -n_t/'^), (10) 

where j'^ is the Lorentz factor of the two fluids observed by the local center-of-mass frame 
of the plasma and m = m+ + m_ . Hereafter, a prime is used to denote the variables of the 
center-of-mass frame. Using these variables, we write 

n±f/^ = -(pU^± ^ J'^) . (11) 
m \ e / 



^The relativistic enthalpy includes the rest mass energy. In a case of perfect fluid gas with specific-heat 
ratio F-i-, it is given by h± — ■m±n± +p±/{T± — 1) +p±. 
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We also define the average and difference variables with respect to the enthalpy density as 



h 
Ah 



n 



n 



— H 

nt 

m 

2m_ 



h_ 2171 , 



n 

"8/1 



m 
2m 



m 



mn 



n_ 



m 



h. 



n 



2m^ 
m 

2m 



m 



(12) 
(13) 

(14) 
(15) 



We find the following relations between the variables with respect to the enthalpy density, 



Ah^ 



h - AfiAh, 
Afih - 



Sfj, 



2fi 



Ah, 



(16) 
(17) 



where /x = m+m_/m^ is the normalized reduced mass and A/i = (m+ — m_)/m is the 
normalized mass difference of the positively and negatively charged particles. It is noted 
th at we have a relation, fi = [1 — (A/i)^]/4. Using the above variables, the same calculations 
of iKoidd (120091 ) yield one-fluid equations from the two- fluid equations ([T]) and ([2]), 



1 

ne 



(pun = 

fih^ 



[ne] 



2ne 
Ah^ 



(18) 
(19) 



ne 2 



fiAh^ 



[ne) 



1 



2ne 



V^{App - Ap) 



ne 



ne 



(20) 



When the relative velocity of the two fluids is not so larg e that t he fri ctional force is propor- 
tional to the relative velocity, according to Appendix of iKoidd (120091 ). the frictional 4-force 
density is given by 

R^" = -rine[J'' - p'^il + e)^]. (21) 

Here, the coefficient t] is recognized as the electric resistivity, Pg is the charge density observed 
by the local center-of-mass frame of the two fluids, p'^ = —U^J^, and the thermal energy 
exchange rate from the negatively charged fluid to the positively fluid is given by 







e {p'^ + J,r){Apn^' + np'Je) 
2ep'^ (n + App'J{2e))n^^ 



(22) 
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where 6 is the redistribution coefScient of the thermahzed energy to the positively and 
negatively ch arged fluids with the equipartition principle (0 < 6* < 1) (see Appendix A of 
Koidel J2009h ) and 



n 



e V e 



(23) 



Using an equation derived by Maxwell equations 



and Equation (jlj), we write the equation of motion (Equation (fT9|) ) by 

V^T'^^ = 0, 

where 



(24) 



1 



\ [neyh J ne 4 

(25) 

This equation corresponds to the eq uation of mo t ion in the ideal GRMHD, for exam- 
ple, Equation (A2) in Appendix A of iKoide et al.l (120061 ). The newly additional term in 
the equation of motion is only that of the current momentum density nh^J^J^/lne)"^ and 
[2/iA/i/(ne)](f/^J^ + J'^f/^) in Equation (|25l). 



To check causahty of Equations f|T8l) - fl20|) . ([3]), and (jl]), we derive the dispersion relation 
of the electromagnetic wave in a uniform, unmagnetized plasma from these equations as 



H 



rik\ fV^\'^ 



(26) 



where k is the wave number, u) is the angular frequency of the electromagnetic wave, and H = 
{rine)"^ /[fi{h^ — 4:ii{Ahy / h}] (see Appendix B). This equation is mathematical ly identical t o 
the di s persio n relation of electromagnetic wave in the resistive pair plasma (IKoidd 120081 ). 
Koidel (120081 ) showed that the group velocity of the electromagnetic wave is smaller than the 
speed of light when H < 2, while it can be larger than the light speed unphysically when 
H > 3. Appendix C shows the inequality 



In A 



16tT J ;U3/2(l/2 



crit ' 



(27) 



where A^p is the plasma parameter and ( = nj^n^/{nj^ + ra_)^ is a variable related to charge 
neutrality. When A^p > A"crit/v^, we have H < 2. This confirms that Equations (|T8l) - 
([20]), da]), dH) for the plasma (A'p > 1 and then A'p > A'crit/v^) are causal. We call these 
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causal equations the "generalized GRMHD equations". Especially, Equation fl20l) is called 
the "generalized general relativistic Ohm's law" . 

When we apply the covariant form of generalized GRMHD equations to cold, quasi- 
neutral plasmas {p± ^ m+n+, ^ n _, i.e. Ah <^ h), Equations (fT8|) - (!20|) are identical to 
Equations (140)- (142) of lKoidel J2009h as, 



0, 



1 

ne 



-V. 



h ( U^'W + 



fx 



(28) 
(29) 



ne y 



A/i 



ne 



-^^^(A/xp - Ap) + (u" - - rf[J^' - p'^{l + Q)U'' 

2ne \ ne J 



(30) 



2.2. 3-|-l formalism 



To understand these generalized GRMHD equations intuitively, we derive a 3+1 formal- 
ism of the equations. We assume that off-diagonal spatial elements of the metric g^^^y vanish, 
gij = {i j). Writing non-zero components by 



9oo — —^0' 9ii — hi , 



we have 



ds^ = g^^dx^'dx" = -hldt^ + ^ [hiidx' f - 2h1ujidtdx^ 



i=l 



When we define the lapse function a and shift vector by 



a 



1/2 



/3' 



hiUi 



a 



the line element ds is written by 

3 

ds^ = -a^df + ^(/ijcix* - a(3'dt) 



i=l 



(31) 
(32) 



(33) 



(34) 



The determinant of the matrix with elements g^^, is given hj g = —{ahih2h^Y , and the 
contravariant metric is written explicitly as 



g'' = 9'' = 9'' = 



hihj 



(35) 
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where 5*-' is the Kronecker's 5 symboL 

We introduce a local inertia frame called the "zero-angular-momentum observer (ZAMO) 
frame". Using the coordinates of the frame {t,x^,x'^,x^), the line element is 

ds"^ = -dP + '^{dx'Y = r]^^dx^dx\ (36) 

i 

where 

dt = adt, (37) 
dx' = hidx' - aP'dt. (38) 

This is identical to the Minkowski space-time locally. In the Boyer-Lindquist coordinates, 
when we write any contravariant vector by a^, according to Equations fl371) and (l38l) . the 
contravariant vector in the ZAMO frame, ci^, is given by 



a° = aa°, a' = Ka' - ajS'a^. (39) 



A covariant vector d^ is 



1 „ 1 

ao = + 2_^-^ai, ai = —ai. (40) 

Note that, because the metric is Minkowskian, we have a° = —So and d^ = di. 

The contravariant and covariant components of vectors and tensors measured by the 
ZAMO frame are given by Equations fl5^ and fHUl) . Denoting these components observed 
by the ZAMO frame with hats, we have 

(41) 

[/° 

(42) 

7 

(43) 
(44) 
(45) 



7 








v' 








e + 7P 




pi 










= = —Foi = 



j ■' 



^eijkBk = Fij - ——Fij, (47) 
Pe = J^ = aj\ (48) 



J* = hiJ'-ap'J^. (49) 
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The relationship between the variables measured in the ZAMO frame is similar to that 
of ideal sRMHD but not identical (iKoide et al.lll996l ). Here, we summarize the relations, 



1 



7 



i\2 



P' = h 

e = h 
= pS'<> + h 



Zneh [nej'^h 



neh (ne)'^h 



+ (Ex B)„ 



Zneh [neYh 



(50) 

(51) 
(52) 



(53) 



where = ^2 ^ ^2 ^ ^2^ ^2 = ^2 _^ ^2 ^ ^2^ 

The generalized GRMHD equations except for the Ohm's law ([TS]) . ffT^ . 
are written as, 



I d_ 



iV^pU") = 0, 



d 



-g dx 



d^E^x + duExn + d\Efj_i, 
1 d 



0, 
0, 



and 



(54) 

(55) 
(56) 
(57) 



where we used the following relations, V^a'' = S^a*^ + rj^^a°", FJ^^ = \g^^{—dpg^a + f^^jfl'po- + 
dadfip), rj^cr = 9^ (In y/—g), and F^,^ = —Ey^. With respect to the Ohm's law (!20|) . defining 
an energy-momentum tensor of charge and current (electric energy-momentum tensor) 



ne 



neh^ 



Ah 



(58) 



we have the following form. 



ne 



1 

ne 



-.dy{^K^')+Ei^yK° 



-^V^(A/xp - Ap) + (u" - ^r) - r][J>' - p'^il + e)U^]. (59) 



-lo- 



using the ZAMO variables, we obtain the following set of equations from equations (|5 
and (1591) 



dip 

dt hih2h 



1 ^ d 



ahih2h-. 
hi 



(60) 



d 



ahih2h^ 



(e + 7p) 



1 da 
hi dx^ 



(61) 



de 

at 



d 



ahih2h^^ 
hi 



j j,k j,k 



(62) 



ne dt \ne J ne 



2fih^ 1 da 



+a 



ne hi dx^ 
1 1 d 



1 d 
hih2h3 ^ dx^ 



hi 



ne 



^ a iGijK'^ - Gjik" + (3^^ (GijP' - GjiP^ | + E —^Ji^^^ 



Ap 



2ne a.^^^^^^ - A^) + - ) " ^^^^ " ^^^^ + ^^^^^ 



(63) 



2^d_ 

ne dt 

+E 



pI + - ( A/ip - Ap) 



ne 4 



1 

ne 



1 >^ d 

hih2h3 ^ dx^ 



ahih2h-i fih^ 
hj ne 



{P + 2/5^pt) 



p/i-'- 1 9a 
ne hi dx^ 



-a 



2ne ^-^ hi dx^ 

dBi -hi 



d 



dt hih2h?j ^ dx^ 



Lm 



AAA A^i \ A J ) ^' 



hih2h2,dx^ V 



(64) 

(65) 
(66) 
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E 



1 d f hih2h 



^1^2 ^3 dx^ 



h, 



-E, 



(67) 



dt 



hi 



ijk 



d 



^1^2 ^3 dx^ 



ahk LBfc + ^ ekimP^Ek 



Lm 



(68) 



where = -^^0, and a^,- _ g^. 



^ ^ 'a(3'^). Here, we used formulae about covariant 



derivative of a tensor (lAlip and (]A12p in Appendix A. This form of the general ized GRMHD 
equat ions is called the 3+1 formalism of the generalized GRMHD equations (IThorne et al. 
19861 ). Here, we also defined following variables. 



ne 
ne 



^00 



-iJ' + pJJ' —pj' + -^^^lU' 



2ph^ 



neAh , 



(69) 

2pht- - 2^^^^ + W^"' 

which can be regarded as modified current density and modified charge density, respectively. 



neh^ 
rPe 



-7 



2.3. Equation of state 



To close the generalized GRMHD equations, we need additional two equations. Here 
we use the equation of state for each fluid. 



= H. ( ^) , (71) 

where ifg is the specific enthalpy, which generally depends on the state and species of the par- 
ticles. In the case of single component relativistic fluid in thermal equilibrium, the function 
is given by 

K,{llx) 



H,{x) 



K^il/x) 



{x > 0), 



(72) 



( IChandrasekharl Il938l : ISyngd 119571 ). Here, K2 and K3 are the modified Bessel functions of 
the second kind of order two and three, respectively. In the case of ideal fluid gas with 
specific-heat ratio F, we have 

^ {x> 0), (73) 



Hs(x) = 1 + 



r - 1 



X 



as shown in the footnote [T] in section [2l Using the equations of states, we find the following 
equations with respect to the variables of the generalized GRMHD equations, 

'p + Ap\ m\ ^ /p — Ap\ m?i 



h 



n 



(74) 
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Ah 



2n fim 



p + Ap \ 
2p+ / P+ 



p — Ap\ m_ 

2p- ; 



where 



m± 
m 



e V e / 



.1/2 



(75) 



(76) 



These simultaneous equations are regarded as the equations of states for the generahzed 
GRMHD equations. In the perfect fluid g with the equal specific-heat ratio F = r+ = 

r_, they reduce to 



h 
Ah 



n 



r 



2(r-i) 



2 limn 



1 1 1 

n+ n„ 2(r 



2 
























{( 


7i 





2 2 

mi 



^2 



(77) 
■ (78) 



When the plasma is cold and quasi-charge neutral, that is, p± <^ p±, and |n+ — n_\ <^ n, 
Ah vanishes and Equations (!75l) and (178!) become superfluous. 



3. Distinctive Properties of Plasmas around Rotating Black Holes 

In this section, we investigate distinctive phenomena of plasmas in space-time around 
rotating black holes, ( 3j J Ob Ob Ob ^) = (t, r, 9, (j)). The metrics of the rotating black hole with 
mass M and angular momentum J are given by the Kerr metric. 



hn 



^1 ^2 



I A . 

— smfc', LUi 



UJ2 = 0, 



2rlar , ^ 



where = GM is the gravitational radius, G is the gravitational constant, a = J/Jmax 
is called the black hole rotation parameter, Jmax = GM^ is the angular momentum of a 



maximally rotating black hole, A = — 2rgr + (a^g)^, S 



+ (arg)^ cos^ 6*, and A 



[r^ -|- (arg)^]^ — A(arg)^ sin^ 6*. The lapse function is a = a/ AS/A. The radius of the event 
horizon is = Tg 



which is found by setting a = 0. 



First, we review the difference between the standard resistive GRMHD and ideal GRMHD 
briefly. The standard resistive GRMHD is given by the conservation laws of particle number, 
momentum, and energy. 



V,(pf/^) = 
V 



(80) 
(81) 
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the standard relativistic Ohm's law, 

UT'', = r/[J^ + {U'^J.p^ (82) 

and Maxwell equations ([3]) and (j4]). The difference between the standard resistive GRMHD 
and ideal GRMHD is only concerning to the Ohm's law. The "ideal MHD condition" of the 
ideal GRMHD is given by the Ohm's law with zero resistivity, 

W'F^, = 0. (83) 

Using the ZAMO frame, we write the standard relativistic Ohm's law (Equation 0821) ) by 

^.p:^ = f/-F,, = jE^ + Y,^^jkU'Bk = v[r + {U''l)U% (84) 

and the ideal MHD condition (Equation fl83|) ) by 

U^'F,, = ^E, + Y^ ^iokU'Bk = 0. (85) 

The significance of the resistive term (the right hand side of Equation 0821) ) is evaluated by 
the ratio r]J^ / {U^ Fi^,) ~ 7]J/{JJB) ~ 7]/{UL), where f/, J, B, and L are the typical values 
of |f/|, I J|, \B\, and a scale length of a system. Inverse of the ratio 

= — , (86) 
■q 

is called the magnetic Reynolds number. When the magnetic Reynolds number Sm is much 
larger than unity, we can neglect the resistivity and use the ideal GRMHD. 

It is noted that causality in the standard resistive GRMHD is broken down when an 
electromagnetic wave packet propagates, while that of the ideal GRMHD is hold because no 
electromagnetic wave can propagate in the ideal MHD plasma. The results of the standard 
resistive GRMHD should thus be treated carefully. 

Next, we move on the evaluation of the generalized GRMHD equations. The difference 
between the generalized GRMHD equations and the standard resistive GRMHD equations 
is shown in the terms with ne and O in Equations fl60l) - fl68l) . Then we find the conservation 
equations of particle number fl60|) and Maxwell equations fl65|) - fl68l) are the same. With 
respect to the equation of motion and energy (Equation fl24l) . or fl6Tl) and (|62|) ). the difference 
is found only in the definition of T^'^ (see Equation fl2S]) ). that is, T^'^ of the generalized 
GRMHD contains the additional terms /i/i^»/[(ne)2/i] and [2/iA/i/(ne)]([/^J'^ + J^'U"). 
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The ratios between the additional terms and the leading term W^U" with respect to e + 7p 
f 0°, P' = f and f'^ are 



Pe 



/i J2 



(716)272' (ne)2^^' {neYtfi' 



i7) 



where we assume the term fih^ J'^ / [{neY h] dominates the term [2pAh/{ne)]{U'^J'' + J^W^) 
in the right-hand side of Equation (|25ll . and the characteristic scales of IC/I, \J\, and \B\ are 
written by U, J, B, respectively. Here, we have the relation between the ratios as 

2 



(ne)2 



(ne)2 72 (ne)2 f/2 

Then, we evaluate the contribution of the term pJ^^J'^/{nef' by the two ratios 

y/PPc \fP J 
ne 7 ' ne IJ 

When we assume quasi-neutrality, we evaluate the ratios by 



ne U 



Vp 



U 



«1, 



(88) 

(89) 

(90) 
(91) 



ne 7 n 

where U± is the spatial components of 4- velocity of positively and negatively charged fluids. 

We find the drastic difference between the Ohm's law of the generalized GRMHD equa- 
tions and standard resistive GRMHD equations. To clarify the newly introduced terms in 
the generalized Ohm's law intuitively, we rewrite it using the identity. 



a(3-^Gji + (Tjj 



hj duj^ 



hjhi dxi 



hi dx'^ ^ ' ^ . g^i ' 
and the identity for the Kerr metric, 

P'G.j = 0. 

We obtain the intuitive 3+1 form of the generalized general relativistic Ohm's law, 

1 1 ^ (9 



(92) 



(93) 



ne dt \ ne 



ne /ii/i2^3 



hi 



ne 



1 1 d 

2ne hi dx^ 



^P 9u 



(A/ip -Ap) + {u''- -^r F,, - r][r - /5;(1 + e)u 



ne 



(94) 



1 

ne 



2iih^ 1 da 
ne hi dx 



jPl 



y^hiduj^ph^ p 
^ hi dx^ ne 
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ne dt 



^pl + hAfip-Ap) 
ne 4 



1 

ne 



_J_S^_d_ 

j 



hj ne 



ne hj dx^ ^ 



hk dx^ 



^ E r^l^(^/^^' - ^P) + - - r/[pe - p1(1 + 6)7] 



2ne ^ h-i dx^ 



ne 



(95) 



We list up the difference between the generahzed Ohm's law of the generalized GRMHD and 
the standard Ohm's law of the standard resistive GRMHD equations as follows. 



1. Hall effect: the term 



ne 



A/i 

ne 



•Ei + ^ tijkP Bk j , 

k / 



(96) 



in the left side of Equation ^ 
2. The inertia of the current: the left hand side of Equation (j9 



3. The transport of momentum by current: the first term of the right hand side of Equa- 
tion m 



4. The thermal electromotive force: the term 

1 1 d 



2ne hi dx^ 
in the right hand side of Equation 



(App - Ap), 



(97) 



5. The term with respect to the equipartition of the thermal energy due to the friction 
between two fluids: the term, rjp'^QW, in the right hand side of Equation (p^ . 

6. Gravitational electromotive force: the first term in the last bracket of the right hand 
side of Equation ( l94l) . 

7. Centrifugal electromotive force: the second term in the last bracket of the right hand 
side of Equation flM|) . 

8. Frame-dragging electromotive force: the last term in the last bracket of the right hand 
side of Equation flUl|) . 
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9. Zeroth component of the Ohm's law: Equation fl95l) . It is used for the calculation 
of A/i. In the case of non-relativistic pressure and quasi-neutral plasma, it becomes 
superfluous {h ^ Ah ^ 0). 



Items (l)-(5) are special relativistic effects, which were already reviewed in iKoidd (120091) . 



Items (6), (7), and (8) are general relativistic effects, which were roughly discussed in iKoide 



( 120091 ) ■ Item (6) shows that the charge separation in the gravity causes the electromo- 



tive forc e , whi ch is also found in non-relativistic gravity. This effect was first reported by 



Khannal (Il998l ). Item (7) shows that the centrifugal force on the electric current causes 
the electromotive force, while item (6) is with respect to the gravitation. Item (8) shows 
the frame-dragging effect causes the electromotive force. The normalized shear of frame- 
dragging angular velocity, [fih^ / {ney][hj / {ahi)]{dujj / dx^) , seems to correspond to the factor 
of a{A^/ne)Fij of the Hall effect rather than to the resistivity. It also appears that the shear 
of the frame- dragging angular velocity corresponds to the magnetic field. However, the nature 
of the coefficient {hi/hj)/{du!j/dx^) and Fij is drastically different because Fji = —Fji and 
dujj/dx^ 7^ —duJi/dxK Thus, these correspondences are just apparent and physically mean- 
ingless. Item (9) becomes significant when the thermal energy density becomes comparable 
to the rest mass energy density and Ah becomes significant. In the nonrelativistic pressure, 
quasi-neutral plasma case, we don't need to solve this equation because h ^ Ah ^ 0. 

We evaluate the contribution of the terms of items (l)-(8) shown above by the compar- 
ison with the term U x B, which is the leading term of the ideal MHD condition (183|1 . We 
write the characteristic length and the time of a system by L and r, respectively. Here, we 
use the relativistic plasma frequency cU p = ne / y/uh and cyclotron frequency tOc = eI3n/h of 



the plasma (the derivation is shown in iKoidd (120091 )). For simplicity, we neglect the contri- 



bution of Ah compared to h for the ordering estimation of terms in the equations. The ratios 
of the contributions of the effects (1)- (8) and the significance of item (9) are as follows: 



1. The Hall effect: 



AfiJ ^ \U+-U\ 

ne U U 



2. The inertia of current: when we consider the first term of the right-hand side of Equa- 
tion (jnnD as the leading term of j\ we have the scaling 

J^lJ_ ^ 1 (99) 

where V = U < 1 is the typical absolute value of 3- velocity and 7 is the typical 
value of the Lorentz factor. 
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3. The transport of current momentum: when we consider the first two terms in the 
parentheses of the right-hand side of Equation fl58l) . we can estimate it as, 



/i/i 1 UJ 

{nef L aUB ' ' (uJpLy 



«7e^T^^z-n^- (100) 



4. The thermal electromotive force: 

a 1 \Afip — Ap\ ^ 1 \A^p — Ap\ 
2ne L aUB ^ 2ujpLU h 



(101) 



5. Equipartition of thermal energy due to friction: this contribution is evaluated by 0. 
When AyU = (pair plasma case) or ^ = 0, this is negligible. However, when A/i ^ 
and 6' 7^ 0, this becomes dominant compared to the last term of the right hand 
side of Equation fl82l) . {U"Jy)tj^^, when pe ^ (quasi-neutrality). We evaluate the 
significance of the friction thermal energy equipartition effect in the generalized Ohm's 
law ( IMI) in the case of quasi-neutrality {\p'J\ -C ne). In the case of quasi-neutrality, we 
approximately have ri^ n from Equation (123!) and then we obtain the approximation 
of B from Equation (122|) 

e.^^. (102) 

where J' is current density observed by the plasma center-of-mass frame. Using the 
definition of the magnetic Reynolds number Sm (Equation fISB]) ) and Equations fll02p . 
the ratio of the thermal energy equipartition term compared to U x B term is 

(f^F,^ ^ B ^ 2SMne 2^Sm 

In the inequality (the last relation) in Equation fll03p . we considered the non-relativistic 
current J' < nej ^fji and an inequality f/ < 7. 



6. Gravitation electromotive force: 

2/i/l a^fBHTPe 2/i Pc fi-BH 2/i Pe 1 



(ne)2 oJjB U ne ujc U ne Luc 



(104) 



where gBu = \{l/ahi){da/dx^)\ ~ / {ahi)]{da / dx^)\ corresponds to the gravita- 
tional acceleration. Here, we used an estimation, (^bh ~ V'^g ~ ^/L. 



7. Centrifugal electromotive force: 



neaUB uj'^rL {uj^Lf 



(105) 
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where r = h^hildh^/dx^) ^ L corresponds to the curvature radius of the x^- 
coordinate hue. 



8. Frame- dragging electromotive force: 



(ne)2 hi dx^ aUB lu^LV ' 



(106) 



where uj^^ = {{h^/ {ahi))duj3/dx^\ ~ \{hj/hi){dujj/dx^)\ is the typical differential an- 
gular velocity of the frame-dragging ZAMO frame. 

9. The zeroth component of the Ohm's law: in the case of Ah <ti h, this equation is 
superfluous because we don't need to calculate Ah. This approximation is available in 
the cases of nonrelativistic pressure and quasi-neutral plasmas. 

Equations (l^ - fllOll) . fll04p - fll06l) show that when the microscopic spatial and temporal 
scales of plasmas are much smaller than the global spatial and time scales, i.e., uj~^,uj~^ ^ 
L, r, cjgjj ,gBB~^, these effects are negligible except for items (1), (5), and (9), because 
\Afip — Ap\ < h and q-qb and uj^^l finite over the whole region around the black hole 
including the event horizon. However, it is noted that when we neglect the current inertia 
effect (item (2)) in the limit case of the small microscopic spatial/temporal scales with finite 
resistivity {L,t ^ u~^, rj > 0), causahty is broken because if 3> 1, while the term of the 
current inertia is negligible compared to the J x B term. When the relative velocity of the 
two fluids is much smaller than the bulk velocity of the plasma, the Hall effect (item (1)) also 
becomes negligible. With respect to the contribution of the equipartition of the thermalized 
energy by the friction (item (5)), it is dominant in the quasi-neutral case (G ^ 1, that is, 
^A/x| J/(ne)p ^ 2p'^/ (ne)) compared to the last term of the right hand side of the standard 
Ohm's law fl82|) . When 9 <^ 1, A/i <^ 1, or ^ it is negligible compared to the 

U X B term. On the other hand, when 6^ ~ 1, A/i ~ 1, ^ and the current is 

sub-relativistic, we cannot neglect the contribution of equipartition of thermalized energy. 
The zeroth component of the Ohm's law (item (9)) is negligible in the case of nonrelativistic 
pressure and quasi-neutral plasmas, while in the case of relativistic pressure around black 
holes we have to solve the zeroth component of the Ohm's law (Equation (195!) ) to determine 
the difference of enthalpy. Ah. 

Then, we conclude that when (i) the magnetic Reynolds number is large enough; Sy^ ^ 
7 / ^/i, (ii) the internal energy (thermal energy and kinetic energy of relative motion of the 
two fluids) of the plasma is relatively small compared to the rest mass energy, and (iii) the 
validity conditions of the generalized GRMHD are satisfled, the ideal GRMHD equations can 
be used for the global dynamics of the plasma around the black hole. If we can not neglect 
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the small spatial or temporal scale phenomena, we have to use the generalized GRMHD. 
For example, through a magnetic reconnection, the small spatial scale event i nfluences the 
global dynamics of the plasma and the black hole rotation (IKoide fc Araill2008l ). Recently, a 
number of numerical simulations using the ideal GRMHD have been performed and showed 
important and interesting physics of the plasmas around black holes, especially, with respect 
to extraction of black hole rotatiorial energy arid relativistic jet forma tion by the magnetic 
field near a black hole (IKoide et al.ll2002l : lKoiddl2003l : iMcKinneyl l2006l ) . In spite of the ideal 
GRMHD calculations, these results showed significant formation of magnetic islands. The 
formation of the magnetic islands requires the magnetic reconnection, which is the result of 
resistivity, and thus the formation of such islands is artificial. However, they suggest that 
the magnetic configuration for the magnetic reconnection, where the anti-parallel magnetic 
fields are almost touched, is easily formed in the black hole magnetospheres. In such a case, a 
significant number of magnetic reconnections will be caused. To perform proper simulations 
with magnetic reconnection, we have to employ the generalized GRMHD, at least, in the 
reconnection region. 



Discussion 



We derived the generalized GRMHD equations from the general relativistic two-fiuid 
equations, and showed the difference among the generalized GRMHD equations, the standard 
resistive GRMHD equations, and the ideal GRMHD equations in section [31 In Table 1, we 
summarize the conditions of validity of each set of equations and the effects neglected by 
the conditions. It is noted that the conditions of sRMHD are given by those with the 
limitation, uj^^, Qbu — ^ 0. Table 1 indicates that the conditions of the standard resistive 
GRMHD equations are too severe to be generally applied to the global dynamics of any 
plasma around black holes. When we consider MHD global phenomena, we can naturally 
assume that <^ L,t, (s'bh)"^, {^bu)~^ where L and r are the characteristic length 

and time of the phenomena, respectively. Within this assumption, we understand that we 
can't neglect charge inertia and current momentum in equations of motion and energy. Hall 
effect, frictionally thermalized energy equipartition effect, gravitational electromotive force, 
and zeroth component of the Ohm's law. Equations (fT8!) - (l2Tl) . ( 163!) . and (16^ with the 
assumption yield a set of equations. 



hUf'U'' 



0, 



+ 



(ne 



Ah 



2ne 



(A/ip - Ap)+ iU 



2ne 
A/i 



(107) 
(108) 



ne 



ju ^ ^f^^j^u + ^ [jM _ + 0)f/M] ^ (109) 
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where the tensor (^j^ with respect to the gravitational electromotive force is defined as 
Co = [2^ihyinef]il/K){d\oga/dx'), C°, = [^ihyinef]il/h,)id\oga/dx-^), Cj = C% = in 
the ZAMO frame. We call this set of equations "modified resistive GRMHD" equations. 

We investigate the conditions of vahdity of the modified resistive GRMHD equations in 
several astrophysical situations. We estimate the critical scale length and time of plasmas 

defined by Lent = max(ccjp ^, ccj^^) and Tcrit = max(a;p ^, u;~^) to clarify the validity of the 
modified resistive GRMHD equations for plasmas around black holes of a GRB, X-ray black 
hole binary, Sgr A*, and AGN (Table 2). When L > L^rit and r > Tcru, the modified resistive 
GRMHD equations arc applicable for the plasma phenomena around the black holes. Here, 
we assume the charge quasi-neutrality of the plasmas and neglect the pressure compared to 
the rest mass energy density to estimate the plasma frequency ujp and cyclotron frequency 
cUc- Using SI units, we have 



-P = ,/^ = J:r^. = (110) 



For the electron-ion plasma, they yield 



^r = — , (111) 

where rrii is the proton mass and rric is the electron mass. We have the relations between the 
frequencies of the electron-ion plasma, a;p^', and those of the pair plasma, w^^^'^, ijJc^^^, 

, .pair _ / 1 c-i ^ e-i pair _ '"i , e-i ^ e-i /t i on 

Cp -a/— '^p >^p , -^^c >^c • (112) 



Using cgs units, we have 



P 



^ 1/2 

P 



4.3 X 10^^ [s-^], (113) 
Vlgcm V 

9.4 X 10^ f-^") [s-^]. (114) 



.IG. 

We also estimate the characteristic scales of the space-time around the black holes by 

5bh = —r-^ ~ — ' 115) 
/is duj'i ac 

c^BH = —r^-[ ' ' (116) 

ahi ox^ rs 

where rs is the Schwarzschild radius, rs = 2rg. We require the variables of <^c~'' 
and rs to evaluate the validity of the modified resistive GRMHD equations. It is noted 
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in) 



rs/(ac) > rs/c and then 



is satisfied if g^l^ ^ 



c/to'p~\ c/uj^~\ For estimation of the variables of plasmas around these black holes, we have 
to give the black hole mass Mbh, plasma density p, and magnetic field B of the plasma. We 
have no direct observation yet of them around any black hole, while a number of estimations 
of them have been done indirectly through theoretical models. Here, we employed the data 
set of Mbh, accretion rate M, p, temperatu re T, and B in t he accretion disks of a GRB, 
X-ray binary, Sgr A* , and AGN estimated by iMcKinneyl (120041 ) (references therein) . Table 2 
shows that the spatial and temporal critical scales, Lcrit and Tcrit, are much smaller than the 
black holes radii (rs) and their light transit times (rs/c) in all of the astrophysical black hole 
objects given here. As an example of the minimum of the characteristic scales of phenomena 
of plasmas [L and r) around the black holes, we consider the minimum scales of magnetic 
reconnection in the plasmas. The minimum scales of the magnetic reconnection are roughly 
estimated by the thickness of the current sheet Lqs, which is calculated by the minimum 
scale of the magnetorotational instability (MRI), 



^cs 



A 



MRI 



4 



2VA 

3 n ' 



where Q is the angular velocity of the disk and va is the A l fven v elocity; Q ^ a/ GM-qr/t^ > 
c/rs and va = a/ B'^/ pop (see Chapter 8 in IShibata et al.l ( 19991 )) and the light transit time 
of the current sheet, Lcs/c. The values of the spatial and temporal scales, Lqs and Lcs/c are 
much larger than the critical variables. Lent and Tcrit, respectively. It suggests the vahdity of 
the modified resistive GRMHD equations even in the phenomena in the reconnection regions 
around the black holes. It is noted that causality is broken in the modified and standard 
resistive GRMHD because the current inertia is neglected artificially {H — > oo), while it is 
kept in the generalized and ideal GRMHD where the current inertia is taken into account 
properly {H < 2) and no electromagnetic wave propagates. However, when we carefully 
generate no electromagnetic wave, we would be able to avoid the problem of causality even 
with the modified and standard GRMHD calculations. 

The generalized GRMHD equations are mathematically identical to the general rela- 
tivistic two-fluid equations, since we used no additional conditi on to der i ve th e generalized 
GRMHD equations. To study the black hole magnetospheres, iKhannal (119981 ) first formu- 
lated the general relativistic two-fluid approximation in the Kerr metric. More generalized 
equati o ns de rived from the general relativistic Vlasov-Boltzmann equations were shown by 
Meier! (120041 ). Here , we c o mpar e the generali z ed GR MHD equations presented here with 
the ones derived by iMeierl ( 12004J ) and iKhannal ( 1l998l ) . The generalized G RMHD equat ions 
derived here are almost intrinsically identical to the equations derived by iMeierl (120041 ) , as 
shown below. For cold (intrinsically non-relativist ic) plasraas, th eir 3+1 formalism reduces 
to the equations quite similar to those derived by iKhannal (Il998l ). while we don't need the 
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condition ''y'^ = 1 whi ch was used b y iKhannal (Il998l ). The basic equation s and approximation 
are different between iMeierl ( 120041 ) and ours. For example. iMeierl (12004 ) employed the Boltz- 
mann equation as a basic equation, while we used the two-ffuid equations. Furthermore, we 
employed rather special definitions of the average variables, especially for enthalpy density, 
h = n^[{h^/n'^) + {h-/n'i_)]. In spite of the diff erence o f the approaches, the effects suggested 
by generalized GRMHD equations derived by IMeierl (120041 ) and us are identical except for 
the Ohm's law. For comparison, we regard the term with respect to the 4-velocity "L /" and 
the in ternal (kinetic) energy flux density "i?" in equations of motion and energy of iMeier 
( 120041 ) (Equation (46) in his paper). 



U ®H + H ®U, 



(117) 



as the terms of Equation (fT9|) of the present paper. 



me 



2ne 



:ii8) 



Hereafter, we also use the notation that IMeierl (l2004r) use d (definitions therein). With respect 
to the Ohm's law, the equations derived by Meiei " |2004 ) and us are similar but not identical 
as follows. We cannot find the left-hand side term of Equation (l20j) . 



[ne] 



(119) 



in the Ohm's law of lMeied (120041 ) (his Equations (57) and (58)), and the factor of the resistive 
term of the right-hand side. 



(120) 



vanishes in the Ohm's law of IMeierl (120041 ) (his Equations (57) and (58)) (see also Equation 
(|2T|) in our paper). The term —p'JJ^ in Equation (11201) is required physically, because the 
Lorentz transformation of the simple Ohm's law in the rest frame, E' = rjj' yields the term. 
Equation (11201) also includes the clear difference between two sets of equations related to 



the effect of the thermalized energy exchange between the positively charged fluid and the 
negative one (the term of G in Equation (I2T!) for the Ohm's law). Furtherr nore, the 4- v ector 
related to current density "j'/Z" of the left-hand side of the Ohm's law of IMeierl (120041 ) (his 
Equations ( 57) and (58)) would correspond to [ph} / {nef']J^ and the current density of the 
Hall term in IMeierl (120041 ) is ^^J — p'^U" , while in the present paper, the current density is J^. 
With respect to the Hall term, we should not subtract the current density due to motion of 
net charge from the total current density, J^, because the net charge moving in the magnetic 
field c auses t he Ha ll effect. With respect to the above two points, the generalized Ohm's 
law of lMeierl (120041 ) should be corrected. When we apply our generalized GRMHD equations 
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in the 3+1 formalism to t he cold plasma, we find the excellent correspondence with the 
equations of iKhannal (Il998l ). except for two points in the Ohm's law. The differences in the 
Ohm's law are only concerning with the term of i n Equation (I2TD and the notation of 
"Pe^"^' current inertia term of the Ohm's law of iKhannal (119981 ) (Equation (65) in his 

paper) and J^* in the present paper (our Equation (l63il ). The current inertia in the Ohm's 
law should not be proportional to the net charge density, Pg, because the current inertia is 
significant even in the cas e of finit e net c urrent with zero net charge. With this point, the 
generalized Ohm's law of IKhannal (119981 ) should be revised. Furthermore, in general, we 
have to consider the contribution of the equipartition effect of thermalized energy due to 
friction, which is des cribed by the t erm propor t ional to 0. Most of the differences of the 



equations derived by IKhannal (119981 ) and 



Meierl (120041 ) from ours come from the differences 



of the definitions of the variables except for the points clarified above. 



We discuss the frame-dragging electromotive force in a similar way as IKhannal (119981 ) . 
When we define antisymmetric and symmetric 3x3 tensors as. 



Hi 



Hji), Hfj 



-{Hij + Hji), 



Hij 



h} hj duj 
ane hi dxi ' 



121) 



we can write the spatial components of the Ohm's law (1941) by 



ne dt \ ne 



+a 



1 1 ^ d 

le /i 1/12/^3 dx^ 

1 1 5 _ . , Afi 



ah,h2hs^^,^ ^^^lihj^i^ 



hi 



ne 



ne hi dx^ 



+ - 



ne 



{Af,p - Ap) + (u^ - F,, - p;(l + 



Because we can regard that corresponds to J"''-', we find [p/(ne)]iJ/^ corresponds to 
[A/x/(ne)]Fjj in the Hall term and \[il{ne)\Hfj corresponds to the resistivity tensor r]bij in 
the resistive term formally. Then, it seems that H^- plays a role of the magnetic field related 
to the Hall effect, and Hfj is related to the resistivity tensor. Furthermore, if we recognize the 
correspondence between A/xFij = A[i{diAj—djAi) and ^H^ = [nh^ / {2ane)][{hj / hi){dujj / dx'') 
{hi/hj){duji/dx^)], we may image the correspondence between AfxAj and [fih^ /{2ane)]hjUj. 
However, this correspondence is just formal appearance and has no real physical meaning, 
because these terms of Hfj and Hfj are not independent and should be considered simultane- 
ously, unlike the Hall and resistive terms. Furthermore, the diagonal elements of Hfj vanish 
and can't cause the resistive phenomena like the magnetic reconnection, which changes the 
topology of the magnetic field configuration and infiuences the whole system drastically. 
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The gravitational electromotive force, 



= ((\^,CVX(\^) = ||v(ln„)pti, (123) 
may cause the magnetic reconnection, while frame dragging and cent rifugal el e ctrom otive 



forces never change the topology of the magnetic field configuration. iKhannal (119981 ) dis- 
cussed this electromotive force but didn't remark the possibility of the magnetic reconnection 
due to the electromotive force. As shown in Fig. (H let us consider a situation where the 
anti-parallel azimuthal magnetic field exists beside the equatorial plane. The current sheet 
is thin and localized near the equatorial plane and the current is directed radially. When 
the net electric charge is distributed at the equatorial current sheet locally, the local radial 
electric field is induced by the gravitational electromotive force E = E^^^. When the direc- 
tion of the gravitational electromotive force -Egrv is the same as that of the current density 
J of the current sheet, we can define the positive effective resistivity rjg^, which satisfies 
Egrv = VgrvJ- This effective resistivity can cause the magnetic reconnection. The sign of 
?7grv depends on the charge separation pi and the directions of current and gravity. This 
mechanism shows that the charge can cause the magnetic reconnection in the strong gravity. 
It is noted that when the plasma containing magnetic fields falls freely into the black hole, 
this effect disappears because there is no effective gravity on the plasma and magnetic field. 
The details will be investigated in near future. 

The magnetic reconnection is expected to happen frequently in the black hole magne- 



tospheres as suggested by a number of long-term ideal GRMHD simulations (IKoide et al. 



2000. 


2006; 


McKinnev 


2006) 



it is not sure whether it makes whole numerical results with ideal GRMHD fatal, because 
the magnetic reconnection region, where the resistivity is significant, is very small compared 
to the global scale, and other regions are able to be treated by the ideal GRMHD. In such 
a case, causality would not be a so serious problem in the small region, and we would be 
able to use the ideal GRMHD for such situations. This speculation should be verified more 
carefully in near future. 

The numerical simulation of the generalized GRMHD will be interesting. The equa- 
tions of state (1711) and (175|l provide closures to the generalized GRMHD equations (I60|) - (l68|) . 
Then, the numerical calculation is primarily possible, although it becomes drastically dif- 
ficult compared to the ideal GRMHD simulations. This is because we have to treat the 
displacement current dE/dt in the Ampere's law and the inertia of the current density 
[1/ {ne)]d[{fj.h^ / {ne)}J]/ dt in the generalized Ohm's law explicitly. In the ideal GRMHD 
calculations, the former is taken into account implicitly and the latter can be neglected 
properly. Furthermore, we have to consider the zeroth component of the Ohm's law to cal- 
culate the enthalpy density difference Ah of relativistically hot plasmas around the black 
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hole. Thus, appropriate simphfication of the generahzed GRMHD equations, especially of 
the simplified Ohm's law, is required for numerical study. The modified resistive GRMHD 
may be one of the candidates. 

In this last paragraph of this section, we investigate the possibility of the stationary 
Ohm's law in a uniform plasma 

(^u' - = r][r - p;(i + e)f/'^], (124) 

to consider the simplification of the generalized Ohm's law (Equations fl20l) and fl2T]) ) for 
numerical simulations. The spatial components of the simplified Ohm's law fll24p are 

- E+(u- ^j] xB = r^[J- p;(l + Q)U]. (125) 

\ ne J \ ne J 

The temporal component of Equation (11241) gives 

U - ^j] ■E = r^iU-J'~ p'^jB). (126) 



ne 

Here, we used the relation J° = •yp'^ + U ■ J'. When we observe the quantities in the plasma 
center-of-mass frame. Equation fll26p yields 

^{J'-E')=vp'^e. (127) 

ne 

Using Equation (I125p in the plasma rest frame, (l — ^' ^ ^ B' = riJ', we obtain 

r? J = 1 p„ 



-^\J'\'=ll-^p'\r^p'^Q. (128) 

ne \ ne J 

Assuming charge quasi-neutrality and non-relativistic collision between the two fluids, we 
substitute the approximation of given by Equation (11021) into Equation (I128p . then we 
find 

—ri\J\ = VPeir^-^ = ir^v\J \ ■ 
ne 2ep\, n 2ne 

This yields 6 = 2 when Ap 7^ and Ohm heating is finite {7]\J'\'^ > 0), while 6 should be 

less than unity. This strange result even for the uniform plasma may suggest the stationary 

Ohm's law contains a self-contradiction. This may come from the unbalance between thermal 

energy gain of the two fluids; otherwise, the non-dimensional variable may have to be 

determined to keep its consistency. This should be clarified in near future. 



I am grateful to Mika Koide, Takahiro Kudoh, Masaaki Takahashi, David L. Meier, and 
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the Science Research Fund of the Japanese Ministry of Education, Culture, Sports, Science, 
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Table 1. Necessary and sufficient conditions for approximation of the generalized 
GRMHD. the standard resistive GRMHD. and the ideal GRMHD approaches. 



approach 


conditions 


effects neglected 


generalized 
GRMHD 


{H<2) 

(definition of the plasma, 1) 




standard 

resistive 

GRMHD 




• current inertia in the Ohm's law 

• current momentum transport 
in the Ohm's law 

• thermal electromotive force 

• centrifugal electromotive force 




1 < u 

< n 


• charge inertia and current momentum 
in equations of motion and energy 




<[/ 


• the Hall effect 


^A/i ^ « 


A 

ne 


• frictionally thermalized energy 
equipartition*^ 


uj^ ^ < (gBu) \ l^pe< net! 


• gravitational electromotive force 




• frame-dragging electromotive force 


Ap ^ P, P'e^ ne {Ah < h) 


• zeroth component of the Ohm's law 


iVp > A^crit/V2 

(condition of generafized GRMHD) 




ideal 
GRMHD 


and all the other conditions in the table 


• resistivity 



^As far as we consider the plasma (A^p ^1), this condition is satisfied well because the plasma 
parameter is much larger than unity. This means that no additional condition of the relativistic two- 
fluid model is required. The conditions of the two-fluid model are given by z/~^, ^ r, Ac, Aj ^ L, 
where u^c and are the electron-electron and ion-ion collision frequencies and Ac and Ai are the 
Debye lengths of electron and ion. For the electron-positron plasma, the ion-ion collision frequency 
and ion Debye length are replaced by the positron-positron collision frequency and positron Debye 
length, respectively. 

^If it is not satisfied, we can't neglect the frictionally thermalized energy equipartition term with 

e. 
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Table 2. Microscopic variables {uj^~\ ■ ■ ■), critical scales (Lent and Tcrit), and 
characteristic scales of phenomena of plasmas around black holes {g^^, ^) 
examine validity of the modified resistive GRMHD equations. 





GRB 


BH X-ray binary 


bupermassive 
BH in Galaxy 


AGN 




GRB030329 


LMC X-3 


Sgr A* 


M87 


Mbh [Mq]- 


3 


10 


2.6 X 10^ 


3 X 109 


Ma 


O.IMq s-1 


IO-^Mq yr-1 


lO-^M© yr-i 


lO-^M© yr-i 


p [gcm~^]'^ 


1.6 X 10^° 


0.0072 


1.5 X 10-6 


1.3 X 10-8 


T [KY 


1.2 X 10^° 


1.8 X 10^ 


3.5 X 10^ 


7.1 X 10^ 


B [Gf 


2.7 X 10^^ 


1.1 X 10^ 


1.4 X 10^ 


3.7 




5.4 X 10^1 


3.6 X 1015 


5.3 X 1013 


4.9 X 1012 




2.5 X lO^s 


1.0 X 101° 


1.3 X 106 


3.5 X 10^ 


[s] 


1.9 X 10-22 


2.8 X 10-16 


1.9 X 10-1^ 


2.0 X 10-13 


!/<-; [s] 


4.0 X 10-19 


1.0 X 10-1° 


7.7 X 10-9 


2.9 X 10-5 


c/Wp"' [m] 


5.6 X 10-1^ 


8.3 X 10-« 


5.7 X 10-6 


6.1 X 10-5 


c/u)l~^ [m] 


1.2 X 10-1° 


0.030 


230 


8.6 X 10^ 


Lcvit [m] 


5.6 X 10-1° 


0.030 


230 


8600 


''"crit [s] 


4.0 X 10-19 


1.0 X 10-1° 


7.7 X 10-9 


2.9 X 10-5 


^BH ~ [m] 


9.0 X 10^ 


3.0 X 10^ 


7.8 X 109 


9.0 X lOii 


^bh/c~ '^s/c [s] 


3.0 X 10-5 


1.0 X 10-^ 


26 


3000 


L > Lcs [m]b 


480 


10 


2.3 X 10^ 


7.5 X 105 


r > Lcs/c [s]-^ 


1.6 X 10-6 


3.3 X 10-s 


7.7 X 10-5 


2.5 X 10-3 



'^Data are from iMcKinneyl (120041 ) 

^'The minimum values of L and r are estimated by the thickness (Lcs) and the 
light transit time (Lcs/c) of the current sheet caused by MRI in the disk. 
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A. 3+1 formalism of divergence of tensors 



We derive a 3+1 formalism of the covariant equation 



(Al) 



for the arbitrary symmetric tensor S^'^, which satisfies S^" = S"^^. We have the relations 
between any tensor S^'^ and the corresponding tensor observed by the ZAMO frame S'^'^, 



coo _ cOO 
o — -o , 



Using T^j^ = (9^ (log- 



a hi 
), we obtain 



1 



hihj 



(A2) 



1 d 



-g dx^ 



With the symmetry of S^'^ ^ we have 



(A3) 



Using Equations flA2l) and flA4p . we get 



where 



j j,k j ^ 



S' 



00 



= -y^g^'idwoo), 
I 

= g^'id,goo) + yg''\d,goi-digo,), 



C'kj = 9'"id,gko)+9''id,gkk), 
1 

r 



(A4) 



(A5) 



(A6) 
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Using Equations fl3T]) . fl33|) . and fl35|) of /ij, and a, we have 



B] 
D] 



^kj 



^^^■^0 + 7^ [-^^(^^"/^O + ^^l^.-"/^')] -JItI- + diih.aP^)] 



2A 



2 hi hi 



dih' 



(^j a 1.2 
2ah-i 



-d,hl 



^d,{h,aP') - ^d,hl 



f3' 



2ah 



2-dkh'j. 



(A7) 



Then, we get 



j I- J 



- V T^TT^ f^*^' + p'S^^ + P'S''' + + ahiT'S^" = aiW + p'H''), (A8) 

hihj ox^ V y A V / V / 



j 



9 ^00 



d 



dt^ ' hih2h3 ^ dx^ 
1 da 



(A9) 



Subtracting Equation (]A9P multiphed by /3* from Equation (]A9p . we have 



_1 >^ d 



8 - 

dt hih2hs ^ dx^ 
j 



'-9 



1 d 



(AlO) 
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Substituting Equations flA5l) and flA7l) into Equations flA9l) and flA10|l . we finally obtain 



(9 ^ 

St hih2h 



1 ^ 9 

22 /is ^ dx^ 



hj 



j,k j,k 



d_ 

di 



d 



hih2hs ^ dx^ 



ahih2h3 



+ 



1 da 



hj dx 



00 



(All) 



J2 a [GijS'^ - GjiS" + p^{GijS^' - + ^ a^iS^'^ = aH\ (A12) 



B. Dispersion relation of electromagnetic wave in unmagnetized plasma 

We derive the dispersion relation of the electromagnetic wave in a uniform, unmagne- 
tized plasma with p = p, f/^ = (1, 0,0,0), p = p, h = h, Ah = Ah, Ap = Ap, F^^ = 0, and 
= in the Minkowski space-time. When perturbations due to the electromagnetic wave 
to the uniform variables are written hjp = p — p, p = p — p, h = h — h. Ah = Ah — Ah, 
Ap = Ap- Ap, = (0, il), F^^ = F^^ {E = E, B = B), = J^", we have the linearized 
equations in the 3-vector form, 

^p + pV-C/ = 0, (Bl) 
d ~ AaAh d d j„ ~ 2pAh„ ~ 

-h + -^-p,--p + hV-U + -^V-J = 0, (B2) 
ot ne ot ot ne 

- d ~ 2pAh d ~ „ , 

^ftt d^j^2j^d^ ^ -Lv(Af.p - Ap) + £ - „J,(B4) 



[neY dt ne dt 2ne 

h^ d I d ~ I d h^ ~ Ah ~ 

Pe + 7^i^h+-^^(App~Ap) + -^V-J + ^V-t/ = 0(B5) 



2{ney dt 2ne dt 8pne dt {2neY 2ne 

V-E = Pe, (B6) 

= 0, (B7) 

= -Vx^;, (B8) 

J+-^E = VxS. (B9) 

In this Appendix, we note an equilibrium variable with a bar, and a perturbation with a 
tilde. We also assume the resistivity is uniform and constant, and the perturbation of any 
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variable A is proportional to exp{ik ■ r — iut) = exp{irifj,i,k^x'^), where = {00, k) is the 
constant contravariant vector called the wave number 4- vector. For simplicity, we investigate 
the transverse mode of the electromagnetic wave in an unmagnetized plasma, thus we set 



E, B, U . 

Then, we have the following linearized equations: 



h + 
—A 



p 



ne 



4/i — pe + Ah 

\ne ^ 



ne 



-iuj-^(—J + 2Ahtj\ 

ne \ne J 

—iojB 
J — iujE 

Pe = 0. 



k. 

= 0, 

- P = 0, 

+ A/ip - Ap = 0, 

= 0, 

= E-r^J, 

= —ik X E, 

= ik X 



(BIO) 

(Bll) 
(B12) 

(B13) 
(B14) 

(B15) 

(B16) 
(BIT) 
(B18) 



After some algebraic calculations, we obtain the dispersion relation of the electromagnetic 
wave, 



T] — lUJ- 



{h}h-4:^i{Ahf) 



(fc^ — oj'^) = iuj. 



(B19) 



{ne^h 

This dispersion equation is identical to that of the electromagn etic wave in t h e pair plasma 
when we reg ard A^{h^-A^{Ahf/h) as the enthalpy density h in koidel hm^ . Ixoidel fl2nn8h 
proved that the group velocity of the dispersion relation (Equation (IBlOp ) is less than or 
equal to the light speed if 

{neniY 



H 



< 2. 



(B20) 



C. Relation between H and plasma parameter 

We show the relation between H and the plasma parameter. Using the definition of 
and Ah (Equations (fT3l) and (fT6ll ). we have 

^i[hh^ - A^i{Ahf] _ 4 h+h_ 

h (n+n_)2^- ^^^^ 
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From Equation flB20l) . we obtain 



n 



h^h- 



(C2) 



where we omit the bars on the mean variables. We consider the nonrelativistic situations 
where the classical resistivity, 

p2 In A 

(C3) 



?7 



27r/im 



is valid (IBellanl |2006| ). where is the average relative velocity of the two fluid particles 
and In A is the Coulomb logarithm, which is an order of 10. We use the expression of the 
averaged velocity by 



P+ 



f-- 



1/2 



(C4) 



where f± is the degrees of the freedom of the two fluid particles. Then, we give the expression 
of the resistivity by 



7] 



In A 

2'7r/im 



P- 



171:71: 



7n_n_ 



-3/2 



(C5) 



Here, we define the plasma parameter of the plasma with/without the charge neutrality by 

ATp = ^K^Xl, (C6) 



where Ad is the Debye length (see iBellanl |2006| ) . 



X2 



P+ 



P- 



(C7) 



In general, we define the plasma as a particle ensemble where the Debye cube contains a 
plenty of both positively and negatively charged particles. 



n+Ao > 1, n_AD > 1. 



(a 



Therefore, we use the condition of the plasma A^p ^ 1 as far as we consider the plasma. 
Using the inequality of the arithmetic mean and geometric mean 



P- 



771:71. 



^71^ 



2 n 



7l_ 

— + — 

P+ P- 



> 



> 



771:71: 



-f- 



P- 



m 1 n_ 



P+ P~ 
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we get 



^j^2 < m+n_ + m_n+ 



2tt 
MA 



-1 



< 



TT 



In A 



/2 



N: 



crit' 



(C9) 



because /± > 1, where ( is the variable related to the charge neutrality, ( = rij^n^ / {nj^+ri-Y . 
For the neutral plasma, C, is 1/4 and C, decreases as the break of the charge neutrality gets 
stronger. Then, if A^'p > Ncrn/V^, we confirm H <2. Here, we evaluate A'crit as 



AT.,,. = t^^,-'/'C'^\ 



2%' 



(CIO) 



Eventually, we have the scaling 

'In A 



0.281/i 



-3/4 



10 



(cii: 



84 (¥) (t^ 



/4 



0.80 (^) U 



/4 



1/4 

1/4 

(for electron-ion plasma, u = 5 x 10 ^) 
-1/4 ■ (C12) 

(for pair plasma, fi = 1/4) 



When A'p is so large that A^p > Ncnt/V^, causality of the GRMHD equations (Equations 
(IISD-lEnD, dSD, and (HD) is satisfied. 
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Fig. 1. — A schematic picture of the magnetic reconnection induced by the gravitational 
electromotive force. A pair of flow is induced by the magnetic reconnection. 



